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Abstract 

Finite  time  thermodynamics,  a  subfield  of  irreversible  thermodynamics,  is  employed  to  model  a  cascaded  thermoelectric  generator, 
which  incorporates  all  the  essential  features  of  a  real  heat  engine.  Control  volume  formulation  of  a  cascaded  thermoelectric  element 
was  carried  out  over  a  small  but  finite  temperature  gap  to  comply  with  the  principles  of  irreversible  thermodynamics.  Three  important 
dimensionless  parameters  are  identified  to  designate  poor  Thomson  effect;  low  thermal  conductivity  and  low  electrical  resistivity  of  a 
good  semiconductor  or  semimetal.  For  ideal  values  of  these  parameters  it  has  been  demonstrated  that  exactly  half  of  the  Joulean  heat 
arrives  at  both  hot  and  cold  junction  when  temperature  maximum  passes  through  the  longitudinal  center  of  the  thermoelectric  element. 
In  general  when  Thomson  heat  cannot  be  neglected,  the  maximum  and  the  minimum  permissible  length  of  any  individual  module  of  a 
cascaded  thermoelectric  device  is  predicted  involving  thermoelectric  properties  of  the  materials  and  passing  current.  It  is  observed  that 
maximum  permissible  length  is  devoid  of  dependence  on  applied  temperature  gradient  whereas  the  minimum  allowable  length  is  not. 
When  half  the  Joulean  heat  affects  hot  end  and  half  the  cold  side,  thermal  conductance  inventory  is  allocated  equally  between  the  high 
and  low  temperature  side  for  best  possible  device  performance.  Finally,  it  has  been  argued  that  the  choice  of  constancy  of  total  conduc¬ 
tance  is  not  only  a  natural  constraint  but  also  a  purely  realistic  design  criterion  for  heat  engines  or  refrigerators.  Such  design  prescriptions 
those  lead  to  the  prediction  of  shape  and  structure  in  any  flowing  system  is  reported  in  open  literature  as  constructal  principle. 
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1.  Introduction 

Thermoelectric  generator  is  a  useful  and  environment 
friendly  device  for  direct  energy  conversion.  Especially 
the  capacity  of  Peltier  and  Seebeck  effect  to  dispense  with 
the  moving  parts  in  the  realm  of  energy  transformation 
from  heat  to  electricity  and  vice  versa  is  more  appealing 
in  such  devices.  With  the  advent  of  semiconductor  materi¬ 
als  the  efficiency  of  a  thermoelectric  generator  can  even  be 
an  alternative  for  the  conventional  heat  engines  [1]. 
Another  perspective  of  thermodynamic  modeling  of  a  ther¬ 
moelectric  generator  is  that  it  includes  all  the  crucial 
features  of  a  real  heat  engine  in  a  relatively  simple  way 
where  closed  form  expressions  are  obtained  for  the  power 
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versus  efficiency  characteristics  [2].  Here  each  generic 
source  of  irreversibility  is  identified  and  quantified  in  this 
process  to  draw  a  one  is  to  one  correspondence  between 
the  conventional  heat  engine  and  the  thermoelectric  gener¬ 
ator.  So,  the  mathematical  modeling  of  a  simple  thermo¬ 
electric  generator  can  also  replace  the  elaborate  task  of 
simulating  an  actual  complex  power  plant,  heat  engine  or 
refrigerator. 

Much  effort  has  been  bestowed  in  finite  time  thermody¬ 
namics  (FTT)  [3-8],  to  model  real  heat  engines.  FTT  mod¬ 
eling  of  thermoelectric  generator  presents  a  full-featured 
analysis  of  real  engines.  Since  all  heat  engine  models  aim 
at  providing  a  realistic  margin  for  an  improvement  of 
actual  systems,  an  analysis  based  on  FTT  figure- of-merit 
hints  a  more  practical  assessment  of  a  maximum  attainable 
improvement  in  comparison  to  the  margin  based  on  Car¬ 
not  efficiency.  Thus,  the  FTT  modeling  is  a  worthy  endea¬ 
vor.  However,  it  is  to  be  noted  that  FTT  modeling  does  not 
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Nomenclature 

A 

surface  area  of  the  heat  exchanger 

Ke 

electrical  conductivity 

C 

some  finite  constant,  Eq.  (41a)  or  (41b) 

K{ 

lattice  thermal  conductivity 

CA 

some  finite  constant,  Eq.  (41d) 

2 

a  dimensionless  parameter,  Eq.  (6b)  or  (15) 

Cb 

some  finite  constant,  Eq.  (41f) 

A 

a  dimensionless  parameter,  Eq.  (6a) 

Ck 

some  finite  constant,  Eq.  (41c) 

e 

nondimensional  temperature,  Eq.  (4a) 

Cz 

some  finite  constant,  Eq.  (41e) 

6 * 

nondimensional  temperature  distribution  with¬ 

Fh 

fraction  of  Joulean  heat  affecting  high  tempera- 

out  Thomson  heat  consideration,  Eq.  (19) 

ture  heat  source,  Eq.  (18) 

0* 

nondimensional  temperature  distribution  with 

Fi. 

fraction  of  Joulean  heat  affecting  low  tempera- 

Thomson  heat  consideration,  Eq.  (8) 

ture  sink,  Eq.  (18) 

P 

electrical  resistivity  of  the  material 

I 

total  current 

o 

electrical  conductivity  or  reciprocal  of  electrical 

Jx 

electrical  current  density  vector  along  x-direc- 

resistivity  of  the  material 

tion 

X 

Thomson  coefficient  of  the  material 

K 

thermal  conductance,  Eq.  (25b) 

e 

dimensionless  length  of  the  device  arm,  Eq.  (4b) 

Kh 

thermal  conductance  of  high  temperature  side 

r 

location  of  maximum  temperature  with  Thom¬ 

heat  exchanger 

son  heat  consideration,  Eqs.  (9),  (11a)  and  (lib) 

kl 

thermal  conductance  of  low  temperature  side 

£* 

location  of  maximum  temperature  without 

heat  exchanger 

Thomson  heat  consideration,  Eqs.  (21)  and  (22) 

L 

length  of  the  leg  of  a  thermoelectric  device 

£ 

dimensionless  temperature,  Eq.  (20) 

n 

a  type  of  semiconductor  material 

P 

type  of  semiconductor  material 

Subscripts 

<2h 

heat  transfer  rate  from  the  high  temperature 

e 

quantities  of  electrical  origin 

source 

1 

quantities  of  lattice  thermal  origin 

<?h 

heat  flow  rate  to  the  hot  end,  Eq.  (12) 

H 

quantities  related  to  high  temperature  source 

2*h 

heat  flow  rate  to  the  hot  end,  Eq.  (23) 

J 

quantities  related  to  component  of  Joulean  heat 

fij 

Joulean  heat  transport  rate,  Eq.  (14) 

L 

quantities  related  to  low  temperature  sink 

Qk 

conducted  heat  transport  rate,  Eq.  (13) 

max 

maximum 

<2l 

heat  transfer  rate  from  the  low  temperature  sink 

min 

minimum 

R 

total  electrical  resistance,  Eq.  (25a) 

n 

n-type  material 

T 

temperature  distribution  function,  Eq.  (1) 

P 

p-type  material 

T^c 

low  temperature  level  at  which  the  device  actu¬ 

X 

along  x  coordinate  direction 

ally  receives  the  heat 

* 

quantities  pertaining  to  no  Thomson  heat  con¬ 

Tnc 

high  temperature  level  at  which  the  device  actu¬ 
ally  receives  the  heat 

sideration 

AT 

applied  temperature  gap  across  a  module 

Superscripts 

U 

overall  heat  transfer  coefficient 

* 

quantities  pertaining  to  Thomson  heat  consider¬ 

X 

coordinate  direction 

ation 

z 

figure  of  merit,  Eq.  (39a) 

0 

constancy  of  physical  parameter 

Greek  symbols 

Symbols 

a 

Seebeck  coefficient  of  the  material 

( > 

represents  an  averaged  quantity 

8 

effectiveness  of  heat  exchanging  equipment 

A 

represents  a  difference 

y 

cost  of  unit  conductance 

K 

thermal  conductivity  of  the  material 

stipulate  the  highest  ceiling  for  the  efficiency  rather  only 
dictates  the  lower  bound  of  the  optimal  efficiency  of  a  heat 
engine  performance  affected  by  finite  heat  transfer  rate  irre¬ 
versibility  [9].  In  practice,  heat  engines  can  operate  between 
the  two  extreme  limits:  one  is  the  reversible  or  maximum 
efficiency  operation  and  another  is  the  irreversible  or  max¬ 
imum  power  condition.  However,  in  practical  situation  the 
optimum  design  criterion  is  a  compromise  between  the  effi¬ 


ciency  and  power  output.  In  the  terminologies  of  thermo¬ 
economics  [10],  optimum  operating  point  is  a  trade  off 
between  the  cost  of  fuel  and  the  cost  of  installed  hardware. 

In  a  FTT  model  generally  all  possible  irreversibilities  are 
attributed  only  to  the  heat  transport  process  external  to  the 
engine  and  not  to  the  internal  conversion  of  heat  into 
power  [11].  For  a  FTT  model  of  a  thermoelectric  generator 
internal  irreversibilities  remain  in  series.  Bypass  heat  leak 
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[12]  incorporated  into  the  modeling  is  an  additional  shunt 
among  other  possible  alternatives  [13-15]  that  make  the 
engine  to  operate  irreversibly.  In  the  present  study,  bypass 
heat  leak  is  identified  as  a  major  contribution  to  the  mea¬ 
sure  of  internal  irreversibility.  The  conducting  mechanical 
support,  which  is  the  locus  of  heat  transfer  across  a  finite 
temperature  gap,  is  the  geometrical  path  of  irreversibility 
transport.  The  bypass  heat  leak  phenomenon  retains  all 
the  essential  features  of  irreversibility  of  the  engine  and 
offers  an  elegant  mathematical  perspective  for  engine  mod¬ 
eling.  Though  in  a  thermoelectric  generator  Joulean  heat¬ 
ing  itself  remains  as  an  inherent  source  of  irreversibility, 
bypass  heat  leak  has  normally  higher  orders  of  magnitude 
than  internal  irreversibility  alone  in  the  range  of  optimum 
engine  performance. 

The  architecture  of  optimized  flow-system,  in  general,  is 
a  commonplace  occurrence  in  engineering  and  nature. 
Solutions  of  many  challenges  have  been  unified  under  the 
single  encompassing  physics  based  theory,  the  constructal 
principle  [16],  which  conceives  that  the  geometry  (shape 
and  structure)  is  generated  in  pursuit  of  global  perfor¬ 
mance  subject  to  global  constraints,  in  flow  systems  the 
geometry  of  which  is  free  to  vary.  In  this  context,  a  thermo¬ 
electric  device  can  be  thought  of  a  flowing  system  through 
which  heat  and  current  flows.  Calculation  of  efficiency  of 
thermoelectric  device  is  reported  in  open  literature  [17]. 
There  the  conditions  and  consequences  of  heat  transport 
between  the  heating  and  the  cooling  medium  and  the  junc¬ 
tions  are  not  addressed.  The  objective  of  the  present  inves¬ 
tigation  aims  at  reporting  finite  time  irreversibility  of  heat 
transport  mechanism;  the  distribution  of  Joulean  heat  into 
the  hot  and  cold  space  and  its  consequences  on  optimal 
allocation  of  heat  exchanger  inventory  and  finally  to  pre¬ 
dict  the  geometrical  shape  and  size  of  each  individual  mod¬ 
ule  of  a  cascaded  thermoelectric  device.  As  it  is  pointed  out 
in  a  treatise  [18]  that  different  assumptions  can  lead  to  dif¬ 
ferent  results,  the  role  of  assumptions  in  describing  the 
model  will  also  be  stressed. 

2.  The  physical  model 


Th 

\7 


Oh 


junction  design,  hot  ends  of  the  two  legs  are  both  electri¬ 
cally  and  thermally  connected  through  a  highly  conductive 
material.  Thermodynamically  this  arrangement  is  equiva¬ 
lent  to  that  of  simple  design  in  which  n-  and  p-legs  are 
joined  end  to  end.  The  lateral  surfaces  of  both  the  legs 
are  insulated  electrically  and  thermally  to  prevent  the  con¬ 
tact  from  each  other.  Additionally,  the  cold  ends  of  the  two 
legs  are  either  insulated  only  electrically  or  situated  sepa¬ 
rately  from  each  other. 

In  literature,  thermodynamics  of  irreversible  process  is 
applied  to  a  thermocouple  where  the  legs  may  have  an  arbi¬ 
trary  shape  and  size;  the  composition  may  be  inhomoge¬ 
neous  and  anisotropic  for  the  transport  quantities  and 
the  properties  of  the  materials  are  arbitrary  functions  of 
temperature  field.  Since,  the  maximum  thermal  efficiency 
of  the  device  is  independent  of  the  shape  of  the  leg  of  a 
thermoelectric  element  [19],  in  our  present  study  its  shape 
and  size  is  immaterial.  The  geometry  and  physical  property 
of  the  n-leg  generally  differs  from  those  of  p-leg.  Here  we 
will  cast  the  problem  using  control  volume  approach  along 
with  the  method  of  average  parameters  [17]. 


To  manifest  the  effect  of  electric  current  and  heat 
transfer  irreversibilities  on  the  thermal  efficiency  of  a  ther¬ 
moelectric  power  generator,  we  consider  the  two -leg  assem¬ 
bly  of  the  basic  components  of  a  device  as  is  shown  in 
Fig.  1.  The  hot  junction  is  maintained  at  a  high  tempera¬ 
ture  level  Tnc  and  it  receives  a  net  heat  transfer  rate  Qu. 
Similarly,  the  cold  end  of  the  two-leg  arrangement  is  held 
at  constant  temperature  TLC  such  that  the  net  heat  rejec¬ 
tion  rate  by  it  is  QL.  Ideally,  the  two  legs  are  one-dimen¬ 
sional  conductors  along  which  x  is  directed  from  THC  to 
Tlc.  The  potential  difference  generated  due  to  Seebeck 
effect  causes  to  flow  a  total  current  /  through  the  total  elec¬ 
trical  resistance  R  of  the  elementary  module  of  length  L  of 
a  cascaded  system. 

The  two  legs  n  and  p  are  generally  chosen  to  be  of  dis¬ 
similar  semiconductors  or  semimetals.  In  a  conventional 


3.  Control  volume  formulation  of  a  single  thermoelectric 
element 


We  will  seek  the  temperature  distribution  along  the 
device  leg,  as  it  is  one  of  chief  importance  for  the  evalua¬ 
tion  of  thermal  efficiency  of  the  device.  Under  steady  state 
condition  for  the  divergence  of  the  flux  vector,  total  energy 
remains  constant  along  any  coordinate  direction  of  a 
dimensional  space.  With  reference  to  Fig.  2,  specializing 
along  x-direction  for  each  leg  we  obtain  [20] 


r  d  T 
T  ^Jx  i 
dx 


j2  d 
pJ'~6x 


=  0 


(1) 


where  Jx  is  the  electrical  current  density  vector  along 
x-direction,  T  is  the  temperature  distribution  function,  k 
is  the  thermal  conductivity  of  the  conductor,  p  is  the 
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Fig.  2.  A  cascading  thermoelectric  element  exposed  to  simultaneous 
current  and  heat  flow. 


electrical  resistivity  of  the  conductor,  a  and  %  are  the  See- 
beck  and  the  Thomson  coefficients,  respectively. 

The  solution  of  this  equation  for  temperature  distribu¬ 
tion  demands  a  specification  of  the  dependence  of  a,  k,  p 
and  t  on  x  or  T.  One  viable  approximation  consists  of 
replacing  all  transport  coefficients  by  their  averages  [17]. 
In  this  spirit,  the  first  term  in  Eq.  (1)  drops  out  and  we 
arrive  at  the  equation 

(2) 


The  resulting  equation  takes  the  form 

d  20  dO  ,  „ 

__y!_  +  A_0 

where 


, 

<*} 

and 

( k)AT  ■ 

The  boundary  conditions  transform  into 

0=\  at  £  =  0 

and 


(5) 


(6a) 


(6b) 

(7a) 


0  =  0  at  £  =  1 .  (7b) 

Solution  of  Eq.  (5)  subjected  to  boundary  conditions  (7a) 
and  (7b)  reads  as 


0* 


t+i  ' 

1  —  exp(/t)_ 


exp(^)  + 


exp  (AZ)+± 
exp(/l)  —  1 


(8) 


where  the  symbol  (  )  represents  an  averaged  quantity. 

Before  attempting  to  solve  the  resulting  simplified  equa¬ 
tion,  it  is  to  be  noted  that  the  approximation  method  is 


valid  only  if 

Thc  ~  T LC 

(3a) 

but 

T HC  >  T LC 

(3b) 

such  that  the  temperature  difference  across  the  thermoelec¬ 
tric  element  A  T=  Tuc  —  TL c  >  0.  These  mathematical 
restrictions  are  of  little  practical  interest,  since  for  opera¬ 
tion  of  the  device  at  higher  efficiency,  temperature  differ¬ 
ence  should  be  as  high  as  possible  on  the  whole  of  the 
thermoelectric  device.  On  the  contrary,  for  very  high  tem¬ 
perature,  the  phenomenological  representation  of  irrevers¬ 
ible  process  is  inappropriate.  Hence,  the  assumption  of 
negligible  temperature  gap  is  consistent  with  the  physical 
theory  developed  in  literature  [21].  In  the  real  world  of 
engineering  design,  it  represents  a  cascaded  system,  where 
power  generation  takes  place  discretely  in  successive  stages 
in  series  with  each  other  and  the  power  is  extracted  at  each 
stage.  With  the  increase  of  number  of  modules,  the  temper¬ 
ature  gap  across  any  module  is  reduced  and  the  discrete 
power  generation  mimics  the  continuous  power  production 
from  a  single  module. 

Now,  we  will  nondimensionalize  Eq.  (2)  using 


t-tlc  =r-rLC 
Thc  —  7Tc  AT 


(4a) 


(4b) 


Now  we  would  like  to  locate  the  regime  of  maximum 
temperature.  This  is  an  important  observation  when  we 
mimic  a  thermoelectric  device  with  that  of  heat  engine 
[2].  In  a  finite-time  heat  engine  model  there  is  a  continuous 
variation  of  temperature  from  heat  source  to  heat  sink 
along  the  physical  path  of  energy  transport.  When  both 
the  legs  of  the  thermoelectric  device  is  of  the  same  length, 
the  location  of  maximum  temperature  in  either  of  the  leg 
of  the  thermoelectric  generator  is  obtained  by  setting 
^  =  0,  which  yields 


r  =  ^n 


’exp(zl)  —  1 
A 


}• 


(9) 


Now,  we  would  like  to  prescribe  some  design  conditions 
those  will  cause  the  temperature  maximum  to  pass  through 
the  geometrical  mid-point  of  the  module  of  a  cascaded 
thermoelectric  device.  Each  individual  module  can  be 
thought  of  an  independent  heat  engine  or  one-dimensional 
insulation  system.  For  a  narrow  temperature  gap  across  the 
module,  the  temperature  maximum  passes  through  the 
mid-point  of  the  device  and  it  experiences  a  minimum 
entropy  generation  or  equivalently  maximum  efficiency 
condition  [22].  From  the  definition  (6a)  the  imposition  of 
the  design  criterion  A  — ►  0  leads  to  the  specification  of  max¬ 
imum  permissible  length  of  any  individual  thermoelectric 
module  as 


Anax^.i  (10a) 

\f)  J  X 

Again  the  design  prescription  j  — * ►  0  stipulates  from  the 
definitions  (6a)  and  (6b)  that  the  minimum  permissible 
length  of  the  individual  thermoelectric  module  as 


=  (fl  A T 

(p)  J x  ' 


(10b) 
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A 

Fig.  3.  Location  of  maximum  temperature  in  concurrence  with  Thomson  effect. 


Thus  it  is  to  be  observed  that  the  minimum  length  of 
the  device  arm  depends  on  the  applied  temperature  gap 
whereas  maximum  permissible  length  is  devoid  of  depen¬ 
dence  of  such  imposed  temperature  gradient. 

Eq.  (10b)  as  a  design  criterion  transforms  Eq.  (9)  into 
the  form 


exp(yl)  —  1 
A 


(1  la) 


Eq.  (11a)  is  plotted  in  Fig.  3,  which  shows  that  for  A^O, 
the  relation  between  £  and  A  is  linear.  Expanding  the  right 
side  of  Eq.  (11a)  analytically  around  the  singular  point 
A  =  0  and  then  passing  to  the  limit,  we  have 


lim  f  =  lim  -  In  {  A- 

A^O  0  A 


f[(' 


A2  A 3 

'A+2l+J!+' 


1 


1 

=  2’ 

(lib) 


Thus  Eq.  (lib)  clearly  demonstrates  that  for  A  =  0,  tem¬ 
perature  maximum  passes  through  the  geometric  mid  point 
of  the  conductor  as  the  electrical  current  changes  the  direc¬ 
tion.  As  long  as  Eqs.  (3a)  and  (3b)  are  valid  the  result  ob¬ 
tained  in  Eq.  (lib)  is  physically  realistic.  So,  in  order  to 
construct  a  cascaded  system  the  length  of  the  first  junction 
should  be  half  of  the  total  permissible  length  of  the  assem¬ 
bly  of  the  thermocouples.  The  geometric  mid  point  will  act 
as  a  heat  source  for  the  next  junction  and  so  on. 

Heat  flow  towards  the  hot  end  is  calculated  by  invoking 
Fourier  conduction  law  as 


•  *  dT 

8h  ~  dl 


c=0 


£=0 


Qj 


A  —  exp(yl)  +  1 
_  Aexp(yl)  -  A 


(12) 

where  the  conducted  heat  through  cross-sectional  area  A 
and  of  length  L  is 


Qk  = 


kAAT 

L 


(13) 


and  the  Joulean  heat  source  of  cross-sectional  area  A  and 
length  L  is 


0j  =  PAlAL 


(14) 


such  that  as  A  —>  0 


(15) 


We  calculate  the  ratio 


& 

Qj 


in  the  limit  A  — >  0  in  order  to 


examine  that  what  proportion  of  Joulean  heat  moves  to 
the  hot  end.  On  calculating  the  limit,  using  L’Hospital’s 
theorem,  we  have 


lim 

A^O 


•  * 
eH 

=  lim 

& 

A^O 

A  -  exp(A)  +  1 
A  exp  (A)  -  A 


1 

2* 


(16) 


From  Eq.  (16)  we  observe  that  exactly  half  of  the  Joulean 
heat  proceeds  to  the  hot  end.  The  first  law  of  thermody¬ 
namics  asserts  that  precisely  fifty  percent  of  the  Joulean 
heat  contributes  to  the  cold  junction. 

If  the  electric  current  changes  its  direction  the  Thomson 
heat  also  changes  its  sign.  It  implies  that  if  two  parallel  con¬ 
ductors  of  almost  the  same  geometrical  and  thermoelectri¬ 
cal  attribute  are  placed  in  communication  with  a  single 
reservoir  and  if  the  same  strength  of  current  flows  in  oppo¬ 
site  directions  through  the  two  conductors,  the  Thomson 
heat  so  generated  by  one  conductor  agrees  nearly  with  the 
Thomson  heat  absorbed  by  the  other.  Thus,  the  effect  of 
Thomson  heat  is  almost  nullified.  Hence  the  net  heat  trans¬ 
fer  interaction  with  the  thermal  reservoir  comprises  of 
rejecting  only  two  Joulean  heating  rates  generated  by  the 
two  conductors.  From  Eq.  (6a)  it  can  be  noted  that  the 
Thomson  effect  need  not  be  absent  even  for  a  vanishingly 
small  value  of  the  parameter  A  as  for  any  individual  cascad¬ 
ing  member,  the  passing  current  and  length  of  the  element  is 
small  and  heat  transfer  irreversibility  phenomenon  over¬ 
whelms  Thomson  effect.  Another  mathematical  way  of 
looking  into  the  problem  is  the  imposition  of  the  restrictions 
that  k,  p  and  a  are  constants.  Then  both  (g)  and 
t  =  r(||)  vanishes,  such  that  Eq.  (2)  further  reduces  to 

d2  y 

ko  +  PAl  =  0 


(17) 
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where  the  constant  values  are  indicated  by  the  zero 
subscript. 

At  the  first  sight  it  seems  that  the  foregoing  argument 
waives  the  imposition  of  the  restriction  that  the  tempera¬ 
ture  difference  should  be  small  if  we  ignore  the  origin  of 
Eq.  (17).  It  is  to  be  noted  that  the  constancy  of  these  ther¬ 
mophysical  and  electrical  properties  demand  in  turn  the 
narrow  temperature  range  of  operation  of  the  individual 
element  of  the  device.  Alternatively,  Eq.  (17)  can  be  readily 
obtained  by  applying  first  law  of  thermodynamics  for  a 
control  volume,  where  conduction  of  heat  takes  place  with 
distributed  heat  source,  without  introducing  the  formalism 
of  irreversible  thermodynamics.  Unlike  thermionic  device  a 
cascaded  thermoelectric  element  works  under  narrow  tem¬ 
perature  range  and  hence  we  can  neglect  the  very  effect  of 
radiation  and  convection.  The  distribution  of  temperature 
along  a  thin  conductor  under  the  influence  of  high  current 
involving  radiative  and  conductive  transfer  is  reported  in 
the  literature  [23-26]  from  a  different  perspective.  Contrary 
to  Thomson  heat  consideration  above,  in  this  limiting  case 
we  have  the  liberty  to  formulate  the  boundary  conditions 
as  follows: 


o 

II 

X 

o 

X 

II 

(18a) 

II 

£ 

P 

X 

II 

tX 

(18b) 

and 

Thc  >  ^LC- 

(18c) 

The  absolute  value  of  both  THc  and  TLC  are  but  small. 

Nondimensional  solution  of  Eq.  (17)  subjected  to  the 
boundary  conditions,  Eqs.  (18a)-(18c)  is 

0.  =  (  (19) 


where 


Po  (A-T)  _  0j 
k0AT  Qk ' 


(20) 


Once  again  the  location  of  maximum  temperature  is  ob¬ 
tained  by  setting  ^  =  0  of  Eq.  (19)  and  the  final  result  is 


Passing  to  the  limit  £  — >  oc  in  Eq.  (21) 


we  notice  that  £*  asymptotically  approaches  to  the  finite 
value  f  So  the  numerical  value  of  £  should  be  high  for  tem¬ 
perature  maximum  to  occur  at  the  geometrical  middle  of 
the  conductor.  Thus  even  for  the  idealized  situation  when 
the  thermoelectric  element  behaves  like  a  resistor  under 
the  influence  of  low  current,  the  placement  of  the  second 
junction  begins  at  the  middle  as  if  the  thermoelectric  ele¬ 
ment  were  not  cascaded. 


Heat  flow  inward  the  high  temperature  side  is  given  by 


6  T 

e.„  =  -kaVx 


•  <30* 

n  n 

=  0J 

e=o 

| 

i 

wi 

We  evaluate  the  ratio 


Q* a 
<2j 


in  the  limit  £ 


lim 

2*h 

=  lim 

C— >oo 

0j 

£ — >oo 

1 

2 


1 

2* 


oo  as 


(23) 


(24) 


Eq.  (24)  confirms  that  only  half  of  the  Joulean  heat  goes  to 
the  hot  end.  Energy  balance  states  that  sharply  half  of  the 
Joulean  heat  arrives  at  the  cold  end  side. 


4.  Control  volume  formulation  for  the  complete 
thermoelectric  device 


In  order  to  maintain  a  consistency  with  the  standard 
notation  of  analysis  prevailing  in  the  literature  we  will 
define  the  relationship  between  electrical  resistance  and 
resistivity;  thermal  conductance  and  conductivity  of  the 
thermoelectric  element  introduced  in  the  foregoing  section. 
Electrical  resistance  R  is  related  to  its  counterpart  resistiv¬ 
ity  p  through 


=  pL 
A  ' 


(25a) 


Thermal  conductance  K  is  dependent  on  conductivity  k  as 


(25b) 


For  the  control  volume  formulation  of  the  integrated 
thermoelectric  device  as  shown  in  Fig.  4,  we  employ 
Newton’s  law  of  cooling  [27].  First  law  of  thermodynamics 


Fig.  4.  Schematic  diagram  of  a  cascading  finite-time  thermoelectric  power 
generator  consisting  of  two  differentially  heated  thermoelectric  elements. 
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analysis  neglecting  Thomson  effect  enables  us  to  write 
down  the  following  heat  transport  equations  in  algebraic 
forms  [2].  Finite  time  heat  transfer  rate  to  the  hot  junction 
Qh  is  given  by 

Qh  =  Kh(Th  -  Tw:)  =  y.lTuc  +  K(Thc  ~  TLC)  -  FHI2R 

(26) 

where  K  and  are  the  thermal  conductances  across  the 
reversible  compartment  and  the  hot  junction,  respectively. 
Tu  is  the  temperature  of  the  high  temperature  source  and 
rHC  is  that  of  thermoelectric  element  such  that  THC  <  TH. 
Fraction  of  Joulean  heat  entering  into  the  hot  junction  is 
Fu.  Eq.  (26)  can  be  rearranged  as 

(K  +  Kh  +  oI)Tu c  -  KTU,  -  (KHTH  +  FHI2R)  =  0. 

(26a) 

Similarly,  finite  time  heat  transfer  rate  to  the  cold  junction 
Ql  is  obtained  as 

0L  =  Kl{Tlc  -  rL)  =  yJTu;  +  K ( T HC  -  TLC)  +  FLI2R 

(27) 


=  [K(KHTH  +  KlTl)  +  KHKLTL }  +  yXLTLl 
+  R[Fh(K  +  Kh)  +  FhK]I2  +  ocFlRI 3.  (30) 

Comparing  like  powers  of  /,  we  have  for  the  term  / 2 
K  +  FlKh  =  (K  +  Kh)Fl  +  KFh.  (31) 

Analogously,  rLC  eliminant  of  Eqs.  (26a)  and  (27a)  with  the 
insertion  of  the  expression  for  THC  from  Eq.  (29a)  we  get 

[Wh+Wl)+Wh]  -  0LKHTHI 
+  R[K  +  FhKl}I2  -  aFHRI2 
=  [(K  +  Kh)KHTH+KKLTh]  -  y-KHTuI 
+  R[Fh(K  +  Kl)  +  FlK]I2  -  a FhRI\  (32) 

Equating  similar  powers  of  I  on  both  sides,  we  have  for  / 2 
K  +  FhKl  =  (K  +  Kl)Fh  +  KFl.  (33) 

Any  particular  solution  of  two  identities  (30)  and  (33)  must 
have  general  functional  form  involving  K ,  and  that 
is 

Fh=Fh(^,^h,^l)  (34a) 

and 


where  is  the  thermal  conductance  across  the  cold  junc¬ 
tion.  Tl  is  the  temperature  of  the  low  temperature  sink  and 
Tlc  is  that  of  thermoelectric  component  such  that  FLC  ^ 
Tl.  Fraction  of  Joulean  heat  entering  into  the  cold  junction 
is  Fl. 

Eq.  (27)  can  be  rewritten  as 

KTHc  ~(K  +  Kh~  al) Tlc  +  (. KlTl  +  FlI2R)  =  0.  (27a) 

For  the  Joulean  heat  distribution  it  is  obvious  that 
Fh+Fl  =  1.  (28) 

Now,  the  system  of  Eqs.  (26a),  (27a)  and  (28)  has  four 
variables  Tnc,  TLC,  FL  and  Fn  rendering  single  degree  of 
freedom.  Choosing  FH  to  be  that  degree  of  freedom  we 
solve  for  THc  and  rLC  to  obtain 

=k[(khth  +kltl)  +/2t?]  +  (Kl  -  oI)(KhTh  +FhI2R ) 

HC  K(Kh  +  Kl)  +  (Kh  +  ol) (Kl  -  ol) 

(29a) 

and 

_K[(KhTh+KlTl)+I2R]  +(Kh+oI)(KhTh+FhI2R) 
LC  K{Kh+Kl)  +  (Kh  +  oI){Kl-oI) 

(29b) 

Now,  we  proceed  to  seek  the  possible  set  of  solutions  for 
the  assumed  unknown  variable  FH  or  FL.  Among  many 
other  methodologies  [28],  we  devise  our  own  based  on 
the  symmetry  of  the  problem  as  follows.  Eliminating  FH c 
between  Eqs.  (26a)  and  (27a)  and  providing  an  expression 
for  Flc  from  Eq.  (29b)  we  obtain 

[K{KhTh  +KlTl)  +KhKlTl]  +  a KlTJ 
+  R[K  +  FlKh]I2  +  aFhRI3 


Fl=Fl(K,Kh,Kl).  (34b) 

On  following  symmetry,  we  assume  a  trial  solution  of  the 
form 


F  =  1  (khkl+kkh+kkh\ 

H  2  \KhKl  +  KKh+KKl) 

\(  KhK±  +  2KKk  \ 
2\KHKL+KKH+KKj' 

Employing  Eq.  (35a)  in  Eq.  (28)  we  obtain 
p  _  1  (KhKl  +  KKl  +  KKl  \ 

L  2  \KhKl+KKh+KKl) 

_  1  (  KhKl  +  2KKl  \ 
2\KhKl+KKh+KKJ' 


(35a) 


(35b) 


Substituting  Eqs.  (35a)  and  (35b)  into  the  identity  Eq.  (31) 
we  have  for  both  sides  a  common  expression 


4KKhKl  +  2K2(Kh  +  Kl)  +  K2hKl 
2(KHKL+KKH+KKt_) 


(36a) 


Similarly,  inserting  Eqs.  (35a)  and  (35b)  into  the  other 
identity,  Eq.  (33),  we  obtain  another  common  expression 


C2 


4KKhKl  +  2K2(Kh  +  Kl)  +  KhK[ 
2(KhKl  +  KK„  +  KKl) 


(36b) 


Eqs.  (36a)  and  (36b)  confirms  that  Eqs.  (35a)  and  (35b)  are 
a  set  of  possible  solutions  for  the  identities  (31)  and  (33). 
Further  by  simple  inspection  we  observe  that  Eqs.  (31) 
and  (33)  admit  the  following  numerical  values 


fh  =  o 

and 

Fl  =  1 

(37a) 

Fh  =  1 

and 

Fh  =  0 

(37b) 

Fu  =  \ 

and 

FL=\. 

(37c) 
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Table  1 


Effect  of  some  limiting  combinations  of  conductance  allocation  on  Joulean  heat  distribution 


Range  or  value 

Fh 

Remarks 

K  =  Finite,  KH  =  Finite,  K L  =  Finite 

1 

2 

1 

2 

II 

K  =  Finite,  KH  =  Finite,  K L  =  Finite 

1 

0 

KhKl  =  —2KKl 

K  =  Finite,  KH  =  Finite,  K L  =  Finite 

0 

1 

KuKl  =  -2KKu 

K=  0,  ^H  =  0,  ^L  =  0 

1 

2 

1 

2 

- 

K  =  Finite,  KH  =  0,  KL  =  Finite 

0 

1 

- 

K  =  Finite,  KH  =  Finite,  KL  =  0 

1 

0 

- 

K=  0,  Kh  =  Finite,  K L  =  Finite 

1 

2 

1 

2 

- 

K  =  Infinite,  KH  =  Finite,  KL  =  Finite 

Kn 

kl 

Kh+Kl 

Kh+Kl 

K  =  Finite,  Kh  =  Infinite,  =  Finite 

1  ( 2 K+Kl\ 

1  (  M  \ 

2 \K+Kh ) 

2  \Kh+Kl  ) 

K  =  Finite,  KH  =  Finite,  K L  =  Infinite 

l(  K*  \ 

2  \K+Kh ) 

1  ( 2K+Kf\ 

2  \K+Kh  ) 

- 

K  =  Infinite,  Kn  =  Infinite,  K L  =  Infinite 

0  or  1 

1  or  0 

- 

Now,  we  will  examine  for  what  combinations  of  K ,  KH  and 
K L  these  numeric  values  are  returned  for  the  functional 
relations  (31)  or  (33).  Eqs.  (35a)  and  (35b)  along  with 
Eq.  (37a)  says  that 

KhKl  =  —2KKft.  (38a) 

For  Eq.  (37b)  to  be  tantamount  with  Eqs.  (35a)  and  (35b) 
one  requires  that 

KhKl  =  -2 KKl.  (38b) 

Equivalency  of  Eqs.  (35a)  and  (35b)  with  Eq.  (37c)  de¬ 
mands  that 

Kh=Kl.  (38c) 

Since,  K ,  KH  and  Kh  are  all  nonnegative  quantities,  only 
Eqs.  (37c)  and  (38c)  is  physically  realistic.  Thus  the  neces¬ 
sary  and  sufficient  condition  for  equipartition  of  Joulean 
heat  produced  is  the  equipartition  of  conductance  alloca¬ 
tions  between  high  temperature  and  low  temperature  heat 
sources.  Furthermore,  when  one  of  the  three  conductances 
runs  to  a  very  high  value  leaving  other  two  in  a  moderate 
range,  Joulean  heat  distribution  again  becomes  unequal. 
Table  1  presents  a  compilation  for  the  fraction  of  Joulean 
heat  distribution  for  different  limiting  combinations  of  con¬ 
ductance  allocations. 

5.  Consequences  of  equipartitioned  Joulean  heat 

The  better  thermoelectric  material  for  direct  energy  con¬ 
version  the  higher  should  have  the  value  of  the  dimension¬ 
less  group  zT,  where  T  is  the  average  absolute  temperature 
and  z  is  the  figure- of-merit  of  the  thermoelectric  material. 
This  dimensionless  parameter  for  a  semimetal  or  semicon¬ 
ductor  can  be  expressed  in  general  as  [29] 

zT  =  a2r0  (39a) 

where  a  is  the  electrical  conductivity  and  is  the  reciprocal 
of  electrical  resistivity.  Thermal  conductivity  can  further 
be  treated  as  the  cumulative  effect  of  electrical  conductivity 
Ke  and  lattice  thermal  conductivity  K\  such  that 

(39b) 


For  higher  orders  of  magnitude  of  the  dimensionless  group 
on  the  left  side  of  Eq.  (39a)  in  conjunction  with  Eqs.  (6a) 
and  (6b)  and  the  conditions  A  — >  0  and  £  — >  oo  we  may  stip¬ 
ulate  that 

t  — >  0+,  k  — >  0+,  p  — >  0  +  and  AT  — >  0  +  .  (40) 

Thus,  we  conclude  that  the  present  analysis  is  valid  even  for 
the  finite  temperature  difference  and  the  effect  of  Thomson 
heat  can  be  neglected  for  a  good  quality  thermoelectric 
material. 

Eqs.  (lib)  and  (16)  together  exhibits  an  interesting 
thermodynamic  property  of  a  thermoelectric  element.  For 
temperature  maximum  to  occur  at  the  geometric  center 
of  a  thermoelectric  element,  exactly  half  of  the  Joulean  heat 
approaches  to  the  hot  end  and  the  other  half  to  the  cold 
side.  Similar  observation  is  repeated  through  Eqs.  (22) 
and  (24). 

From  Eqs.  (37c)  and  (38c)  we  learn  that  both  the  hot 
and  cold  junction  experiences  half  of  the  net  Joulean  effect 
produced  and  which  demands  in  turn  equal  conductance 
allocation  on  both  sides.  Eq.  (38c)  also  prompts  the  fact 
that  both  KH  and  KL  are  finite  and  hence  the  following 
proposition  holds 

Kh+Kl  =  C  (41a) 

where  C  is  some  finite  constant.  In  engineering  literature 
conductance  is  denoted  as  a  product  of  overall  heat  trans¬ 
fer  coefficient  U  and  related  surface  area  A.  Thus  Eq.  (41a) 
can  be  rewritten  as 

UHAH  +  UlAl  =  C.  (41b) 

Allocation  of  heat  exchanger  inventory  was  extensively 
investigated  in  connection  with  the  optimization  of  refrig¬ 
eration  and  power  production  both  from  thermodynamic 
[30-32]  and  thermoeconomic  [33]  viewpoint.  But  the  final 
result  of  an  optimization  problem  depends  on  the  nature 
of  imposed  constraint  [18].  Klein  [30]  considered  the  con¬ 
straint  of  the  type 

£h^H  +  =  Cja  (41c) 

where  s  is  the  effectiveness  of  the  heat  exchanging  equip¬ 
ment  and  CK  is  a  constant.  Based  on  the  notion  that  both 


k  =  Ke  +  K\. 
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conductance  and  entropy  generation  have  the  same  dimen¬ 
sion,  Ait-Ali  [31]  conceived  a  condition  of  the  from 


Qh  ,  Ql 
Th  -  rHC  TLC  -  TL 


Ca 


(41d) 


where  Ca  is  some  parametric  constant.  On  the  basis  of 
total  cost  conservation  of  heat  exchanger  installation, 
Antar  and  Zubair  [33]  framed  a  relation  as 

7h^h^h  +  7l^l^l  =  Cz  (41e) 


where  y  is  the  unit  conductance  cost  and  Cz  has  a  fixed 
value.  Bejan  [32]  considered  the  maximization  of  power 
production  of  heat  engines  and  refrigeration  load  in  refrig¬ 
erators  with  two  heat  reservoirs  considering  total  area 
constraint  for  the  heat  exchangers  on  following  the 
equation 

Ah  +  v4l  =  CB  (41f) 

where  A  is  the  surface  area  of  the  heat  exchanger  and  CB  is 
a  constant  due  to  some  resource  constraint.  Treating  Eq. 
(41b)  also  as  a  constraint,  Bejan  concluded  that  in  both 
the  cases  either  quantitative  or  qualitative  equipartition 
of  thermal  conductance  inventory  is  valid.  Similar  results 
have  been  echoed  in  other  works  [31-33].  Thus  we  propose 
that  Eq.  (41b)  is  the  most  natural  constraint  for  such  cate¬ 
gory  of  optimization  problems.  The  idea  of  such  synthetic 
constraint  in  optimization  problems  was  advanced  by  the 
present  authors  [34]. 


6.  Conclusions 

(i)  Creditably,  the  debatable  concept  of  endoreversibility 
[35-37]  in  finite-time  thermodynamics  can  be  mitigated  by 
incorporating  some  irreversibility  factors  to  the  reversible 
compartment  sandwiched  between  two  irreversible  cham¬ 
bers.  Consideration  of  bypass  heat  leak  is  a  compensating 
measure  to  this  direction.  Joulean  heating  present  in  a  ther¬ 
moelectric  generator  itself  is  an  inherent  source  of  irrevers¬ 
ibility  and  tantamount  to  the  frictional  loss  in  a  heat 
engine.  The  discrimination  between  frictional  heat  leak 
and  heat  loss  due  to  finite  rate  heat  transfer  was  first  put 
forward  in  Ref.  [38].  In  the  current  investigation  bypass 
heat  leak  is  identified  as  a  major  contribution  to  the  mea¬ 
sure  of  irreversibility.  The  sufficiency  of  bypass  heat  leak 
consideration  in  engine  modeling  is  an  established  practice 
[2]- 

(ii)  Our  analysis  shows  that  in  a  thermoelectric  genera¬ 
tor  Thomson  effect  may  be  neglected  in  one  limiting  case 
or  may  not  be  negligible  in  another  limiting  situation  even 
for  a  vanishingly  small  value  of  a  certain  nondimensional 
parameter  A.  A  very  high  value  of  another  dimensionless 
factor  £  recognizes  a  better  figure-of-merit  and  the  opera¬ 
tion  of  the  thermoelectric  device  as  cascaded  system  over 
a  small  but  finite  temperature  gap. 

(iii)  Three  parameters  A,  X  and  £  so  identified  are 
responsible  for  temperature  maximum  to  pass  through 


the  geometrical  middle  of  the  one-dimensional  physical 
device.  This  observation  is  in  conformity  with  the  principle 
of  insulation  design  and  the  broader  sense  of  engine  mod¬ 
eling.  The  parallelism  between  the  design  principle  of  heat 
engine,  heat  exchanger  and  refrigerator  to  that  of  insula¬ 
tion  system  was  established  by  a  pioneering  work  of  Bejan 
[22].  For  the  most  efficient  system,  a  stack  of  insulation  is 
cooled  midway.  Similarly  the  calculation  of  mid  point  tem¬ 
perature  is  of  intrinsic  importance  also  in  the  application  of 
thermionic  elements  where  an  interesting  phenomenon 
occurs  in  the  middle  of  the  conductor:  the  temperature 
reaches  extremum,  remains  constant  there  and  Joulean 
heating  and  radiative  heat  transfer  takes  an  equal  share 
of  the  feeded  energy  [23-26]. 

(iv)  Imposing  the  design  prescription  A  — >  0  leads  to  the 
limit  of  maximum  permissible  length  of  any  individual 
module  of  the  cascaded  thermocouple.  Another  design  cri¬ 
terion  j  — >  0  stipulates  the  minimum  permissible  length  of 
such  module.  It  is  to  be  noticed  that  the  maximum  allow¬ 
able  length  depends  explicitly  only  on  the  thermoelectric 
properties  of  the  material  and  the  value  of  the  passing  cur¬ 
rent  whereas  the  minimum  permissible  length  additionally 
depends  on  applied  temperature  gradient. 

(v)  For  the  ideal  values  of  these  three  parameters  A,  A 
and  £  it  is  interesting  to  report  that  exactly  half  of  the  Jou¬ 
lean  heat  flows  into  the  hot  end  and  half  to  the  cold  junc¬ 
tion.  It  is  to  be  noted  that  this  is  not  equivalent  to  state 
that  half  the  Joulean  heat  affects  either  end  [39]. 

(vi)  Control  volume  formulation  for  the  integrated 
device  disseminates  that  exactly  half  of  the  Joulean  heat 
affects  both  the  junction.  This  observation  is  duly  sup¬ 
ported  by  experimental  evidences  [40,41]  in  a  similar  class 
of  thermoelectric  devices. 

(vii)  It  is  conditioned  that  for  finite  bypass  heat  leak 
optimal  conductance  allocation  is  equipartitioned  between 
high  and  low  temperature  side  when  Joulean  heat  affects 
both  the  junction  equally.  This  instance  of  equipartition 
also  conforms  with  the  corollary  [42]  of  constructal  theory. 

References 

[1]  A.F.  Ioffe,  Semiconductor  Thermoelements  and  Thermoelectric 
Cooling,  Infosearch  Ltd,  London,  1957. 

[2]  J.M  Gordon,  Generalized  power  versus  efficiency  characteristics  of 
heat  engines:  the  thermoelectric  generator  as  an  instructive  illustra¬ 
tion,  Am.  J.  Phys.  59  (1991)  551-555. 

[3]  F.L.  Curzon,  B.  Ahlborn,  Efficiency  of  a  Carnot  engine  at  maximum 
power  output,  Am.  J.  Phys.  43  (1975)  22-24. 

[4]  B.  Andresen,  P.  Salamon,  R.S.  Berry,  Thermodynamics  in  finite  time, 
Phys.  Today  (1984)  62-70. 

[5]  A.  De  Vos,  Reflections  on  the  power  delivered  by  endoreversible 
engines,  J.  Phys.  D:  Appl.  Phys.  20  (1987)  232-236. 

[6]  J.M.  Gordon,  Maximum  power  point  characteristics  of  heat  engines 
as  a  general  thermodynamic  problem,  Am.  J.  Phys.  57  (1989)  1 1 36— 
1142. 

[7]  G.  De  Mey,  A.  De  Vos,  On  the  optimum  efficiency  of  endoreversible 
thermodynamic  process,  J.  Phys.  D:  Appl.  Phys.  27  (1994)  736-739. 

[8]  J.  Chen,  The  maximum  power  output  and  maximum  efficiency  of  an 
irreversible  Carnot  heat  engine,  J.  Phys.  D:  Appl.  Phys.  27  (1994) 
1144-1149. 


A.K.  Pramanick,  P.K.  Das  I  International  Journal  of  Heat  and  Mass  Transfer  49  (2006)  1420-1429 


1429 


[9]  Z.  Yan,  L.  Chen,  The  fundamental  optimal  relation  and  the  bounds 
of  power  output  and  efficiency  for  an  endoreversible  Carnot  engine,  J. 
Phys.  A:  Math.  Gen.  28  (1995)  6167-6175. 

[10]  S.  Sieniutycz,  P.  Salamon  (Eds.),  Finite-Time  Thermodynamics  and 
Thermoeconomics,  Taylor  &  Francis,  New  York,  1991. 

[11]  M.H.  Rubin,  Optimal  configuration  of  a  class  of  irreversible  heat 
engines,  Part  I,  Phys.  Rev.  A.  19  (1977)  1272-1276. 

[12]  A.  Bejan,  H.M.  Paynter,  Solved  Problems  in  Thermodynamics,  Dept. 
Mech.  Eng.,  MIT  Press,  Cambridge,  MA,  1976,  Problem  VII-D. 

[13]  1. 1.  Novikov,  The  efficiency  of  atomic  power  stations,  J.  Nucl.  Energy 
II  7  (1957)  125-128. 

[14]  M.M.  El-Wakil,  Nuclear  Power  Engineering,  McGraw  Hill,  New 
York,  1962,  pp.  162-165. 

[15]  P.C.  Lu,  On  optimal  disposal  of  waste  heat,  Energy  5  (1980)  993-998. 

[16]  A.  Bejan,  Shape  and  Structure  from,  Engineering  to  Nature, 
Cambridge  University  Press,  Cambridge,  UK,  2000. 

[17]  B.  Sherman,  R.R.  Heikes,  R.W.  Ure  Jr.,  Calculation  of  efficiency  of 
thermoelectric  device,  J.  Appl.  Phys.  31  (1960)  1-16. 

[18]  A.  De  Vos,  B.  Desoete,  Equipartition  principle  in  finite-time 
thermodynamics,  J.  Non-Equilib.  Thermodyn.  25  (2000)  1-13. 

[19]  A.H.  Boerdijk,  Contribution  to  a  general  theory  of  thermocouples, 
J.  Appl.  Phys.  30  (1959)  1080-1083. 

[20]  T.C.  Harman,  J.M.  Honig,  Thermoelectric  and  Thermomagnetic 
Effects  and  Applications,  McGraw  Hill,  New  York,  1967,  p.  276. 

[21]  S.R.  De  Groot,  Thermodynamics  of  Irreversible  Process,  Interscience 
Pub.  Inc.,  New  York,  1952,  pp.  141-162. 

[22]  A.  Bejan,  Entropy  Generation  through  Heat  and  Fluid  Flow,  John 
Wiley  &  Sons,  New  York,  1982,  pp.  173-177. 

[23]  S.C.  Jain,  Sir  K.S.  Krishnan,  The  distribution  of  temperature  along  a 
thin  rod  electrically  heated  in  vacuo,  Part  I:  Theoretical,  Proc.  Roy. 
Soc.  Lond.  A  222  (1954)  167-180. 

[24]  S.C.  Jain,  Sir  K.S.  Krishnan,  The  distribution  of  temperature  along  a 
thin  rod  electrically  heated  in  vacuo,  Part  II:  Theoretical,  Proc.  Roy. 
Soc.  Lond.  A  225  (1954)  1-7. 

[25]  S.C.  Jain,  Sir  K.S.  Krishnan,  The  distribution  of  temperature  along  a 
thin  rod  electrically  heated  in  vacuo,  Part  III:  Experimental,  Proc. 
Roy.  Soc.  A  225  (1954)  7-18. 

[26]  S.C.  Jain,  Sir  K.S.  Krishnan,  The  distribution  of  temperature  along  a 
thin  rod  electrically  heated  in  vacuo,  Part  IV:  Many  useful  formulae 
verified,  Proc.  Roy.  Soc.  A  225  (1954)  19-32. 


[27]  P.  Salamon,  A.  Nitzan,  Finite  time-optimizations  of  a  Newton’s  law 
Carnot  cycle,  J.  Chem.  Phys.  74  (1981)  3546-3560. 

[28]  K.  Rektorys  (Ed.),  Survey  of  Applicable  Mathematics,  Iliffe  Books 
Ltd,  London,  1969,  pp.  70-75. 

[29]  G.  Min,  D.M.  Rowe,  Thermoelectric  figure-of-merit  barrier  at 
minimum  lattice  thermal  conductivity,  Appl.  Phys.  Lett.  77  (2000) 
860-862. 

[30]  S.A.  Klein,  Design  consideration  for  refrigeration  cycles,  Int.  J. 
Refrig.  15  (1992)  181-185. 

[31]  M.  Ait-Ali,  Maximum  power  and  thermal  efficiency  of  an  irreversible 
power  cycle,  J.  Appl.  Phys.  78  (1995)  4313-4318. 

[32]  A.  Bejan,  Theory  of  heat  transfer-irreversible  power  plants,  Part  I: 
The  optimal  allocation  of  heat  exchange  equipment,  Int.  J.  Heat  Mass 
Transfer  38  (1995)  433-444. 

[33]  M.A.  Antar,  S.M.  Zubair,  Thermoeconomic  considerations  in  the 
optimum  allocation  of  heat  exchanger  inventory  for  a  power  plant, 
Energy  Convers.  Manage.  42  (2001)  1169-1179. 

[34]  A.K.  Pramanick,  P.K.  Das,  Heuristics  as  an  alternative  to  variational 
calculus  for  optimization  of  a  class  of  thermal  insulation  systems,  Int. 
J.  Heat  Mass  Transfer  (2005)  1851-1857. 

[35]  D.P.  Sekulic,  A  fallacious  argument  in  the  finite-time  thermodynam¬ 
ics  concept  of  endoreversibility,  J.  Appl.  Phys.  83  (1998)  4561-4565. 

[36]  B.  Andresen,  Comment  on  “A  fallacious  argument  in  the  finite-time 
thermodynamic  concept  of  endoreversibility”,  J.  Appl.  Phys.  83 
(2001)  6557-6559. 

[37]  D.P.  Sekulic,  Response  to  “Comment  on  ‘A  Fallacious  argument  in 
the  finite-time  thermodynamics  concept  of  endoreversibility’”,  J. 
Appl.  Phys.  90  (2001)  6560-6561. 

[38]  B.  Andresen,  P.  Salamon,  R.S.  Berry,  Thermodynamics  in  finite-time: 
extremals  for  imperfect  heat  engines,  J.  Chem.  Phys.  66  (1977)  1571— 
1577. 

[39]  W.H.  Luke,  Reply  to  experiment  in  thermoelectricity,  Am.  J.  Phys.  28 
(1960)  563. 

[40]  J.H.  Noon,  B.J.  O’Brien,  Sophomore  experiment  in  thermoelectricity, 
Am.  J.  Phys.  26  (1958)  373-375. 

[41]  J.K.  Logon,  J.R.  Clement,  H.R.  Jeffers,  Resistance  minimum  of 
magnesium:  heat  capacity  between  3  °K  and  13  °K,  Phy.  Rev.  105 
(1957) 1435-1437. 

[42]  A.K.  Pramanick,  P.K.  Das,  Note  on  constructal  theory  of  organiza¬ 
tion  in  nature,  Int.  J.  Heat  Mass  Transfer  48  (2005)  1974-1981. 


